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We report on the highly extinguished afterglow of GRB 070306 and the prop- 
erties of the host galaxy. An optical afterglow was not detected at the location of 
the burst, but in near- infrared a doubling in brightness during the first night and 
later power-law decay in the K band provided a clear detection of the afterglow. 
The host galaxy is relatively bright, R ~ 22.8. An optical low resolution spec- 
trum revealed a largely featureless host galaxy continuum with a single emission 
line. Higher resolution follow-up spectroscopy shows this emission to be resolved 
and consisting of two peaks separated by 7 A, suggesting it to be [O n] at a red- 
shift of z — 1.49594 ± 0.00006. The infrared color H — K = 2 directly reveals 
significant reddening. By modeling the optical/X-ray spectral energy distribu- 
tion at t = 1.38 days with an extinguished synchrotron spectrum, we derive 
Ay = 5.5 ± 0.6 mag. This is among the largest values ever measured for a GRB 
afterglow and visual extinctions exceeding unity are rare. The importance of 
early NIR observations is obvious and may soon provide a clearer view into the 
once elusive 'dark bursts'. 

Subject headings: dust, extinction — cosmology: observations — galaxies: red- 
shift and starburst — gamma rays: bursts 



1. Introduction 



Long gamma-ray bursts (GRBs) are believed to be indicators of star formation due 



to their association with short lived massive stars (iGalama et al.l 1 19981 : iHjorth et al.l 12003 



Stanek et al.l 120031 ) . These events should therefore offer grea t prospects for characterizing 
the star formation history of the Universe (jWijers et al.lll998l ). as GRBs can be detected (in 
gamma and X-rays) at all observable redshifts and through large columns of dust and gas. 
One potential objection to this is if GRBs are biased towards a special subset of massive 
stars, i.e. those of low metallicity or in peculiar circumstellar environments. On the other 
hand, if there is no dependency on metallicity then the nature of the h ost galaxies leads to 
the surprising conclusion that most stars are formed in dwarf galaxies (IFynbo et al.l 120071 ). 



1 Based on observations made at the European Southern Observatory, Paranal, Chile under program 
078.D-0416, 177.A-0591; with the Nordic Optical Telescope, operated on the island of La Palma jointly by 
Denmark, Finland, Iceland, Norway and Sweden in the Spanish Observatorio del Roque de los Muchachos 
of the Instituto de Astroffsica de Canarias; with the William Hershel Telescope operated on the island of La 
Palma by the Isaac Newton Group in the Spanish Observatorio del Roque de los Muchachos of the Instituto 
de Astroffsica de Canarias 
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In the pre-Swift era the detection of optical afterglows (OAs) was a requirement for 
sub-arcsec localizations. This was obtained for about 30% of the GRBs due to effects such 
as observability, weather conditions etc. However, in a few cases stringent u pper limits 
on the early (optical) afterglow could be obtained (see e.g. Groot et al. 1998h . providing 
evidence for a sub-class of OA- less bursts dubbed "dark burst" ( see e.g. 



Lamb fc Reichartl I2OO0I : iFvnbo et all l200ll : berger et allbooi lLazzati et al 



Fruchter 1999 



2002af ). Today 



the situation is q uite different as the majority of bursts ( > 90%) are localized with the Swift 
XRT instrument (IGehrels et al.ll2004j ; iBurrows et al.ll2005l ) . IJakobsson et al.l (12004 ) proposed 
an unbiased definition of a dark burst by suggesting that the spectral slope between the 
optical and X-ray is fiox < 0.5. As many as 25% of the Swift bursts fulfill this criterion 
(IFynbo et al.l 120071 ). Throughout the paper we adopt this definition of a dark burst. There 
are several effects that may cause an afterglow to be classified as dark, but among the 
most obvious and p ossibly dominant effects are dust obscuration and high redshift (see also 



Roming et al.l 120061 ) . 



Although the afterglows of dark bursts do not all go by undetected, there is a sizeable 
sub-class of bursts for which we do not have much information and these may consequently 
affect th e luminosity func tion and statistics of GRBs and their host galaxies significantly 



(see e.g. 



Rol et al.ll2007bl . who reported on a dark burst with significant extinc tion). There 



are a lso a few systems which do not resemble the general population of GRBs (jPerley et al. 



2008|). 



In this paper we present the detection of the NIR afterglow of GRB 070306 which is 
characterized as a dark burst. The most likely interpretation is that the afterglow is heavily 
dust obscured and the level of extinction provides a unique opportunity to study the local 
dust properties of the progenitor. The outline of this paper is as follows: in Section [2] we 
describe the Swift X-ray data, and our optical and near- infrared observations of GRB 070306. 
The host and the spectral energy distribution (SED) template fitting is described in Sectional 
the afterglow and the extinction curve fits are discussed in Section HI while the results are 
summarized and discussed in Section [5j 

Throughout the paper we assume a cosmology with H = 70 km s _1 Mpc _1 , fl M = 0.3 
and f^A = 0.7. The afterglow is characterized by a standard temporal and spectral power-law 
on the form, F oc t~ a v~^ . All quoted errors are 1-sigma (68%), unless noted differently. 
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Observations and reductions 



GRB 070306 was discovered by Swift at (to) 16:44:28.0 UT on 2007 March 06 (IPandev et al 
20071 ). The burst had a Tgo = 210 ± 10 s (IBarthelmy et al.l 120071 ). Swift slewed immediately 



and XRT observations began a t t n + 153 s and with U VOT at t + 162 s. An X-ray afterglow 
was identified and reported by lPage &: Pandeyl (120071 ). 



2.1. X-ray 

The Swift-XRT lightcurve shows a typical plateau phase, ending at approximately 35 ks 
after the BAT trigger. To produce the X-ray spectrum at the time of the optical observations, 
the PC events after the plateau phase (35-266 ks) were reduced in a standard way using the 
most recent calibration data and normalised to the count rate at the time of the optical 
observations (to + 1.38 days). A power law with fixed Galactic and variable absorption at 
the redshift of the host (z = 1.4959, see sections [3] and H]) fits the data well (x 2 /dof = 
27.1/31). The best-fit photon index T = 2.2 ± 0.1 (fix = 1-2 ± 0.1); the excess absorption 
iYn = (2.6 ± 0.4) x 10 22 cm -2 for the equivalent hydrogen colum n density (assuming solar 



abundances), in excess of the Galactic value (3 x 10 20 cm -2 , iDickey fc Lockmanl Il990l ) 
corresponding to (0.4±0.04) x 10 22 cm -2 at z — 0. We note that the fit with excess absorption 
at z = is poor. For sub-solar (Z = 0.1Z Q ) metallicity we find 7V~h = (20 ± 3) x 10 22 cm -2 . 
There is no evidence for a change in the spectral parameters between the plateau and post- 
plateau phases: the best fit absorption and spectral slope parameters were consistent to 
within one sigma. 

The optical to X-ray spectral slope can be estimated at Mar 07.00 UT by using the 
specific 3keV flux at this time (0.79 /xJy) and the estimated upper limit on the afterglow 
from the I band NOT image at Mar 06.95 UT (/ > 22 mag(Vega), or F v < 4.1 /iJy). This 
gi yes an upper limit of (3p x < 0.23. This is, hence, a dark burst according to the definition 



of IJakobsson et al.l (120041 ). which requires that Pox < 0.5. 



2.2. Optical and near-infrared 



We carried out early follow up to localize the afterglow in the optical using the 2.5 m 
Nordic Optical Telescope (NOT) and ALFOSC instrument and in the near-infrared using the 
4m William Herschel Telescope (WHT) and LIRIS instrument on Observatorio del Roque 
de los Muchachos. Observing conditions were good, but the M oon was close to being full. 
Imaging and photometry data from the SDSS (jCool et al.ll2007l ) revealed an object close to 
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the detection limit within the error circle. iRol et al.l (j2007al ) reported the detection of a fairly 
bright near-infrared (NIR) after glow detected 3.3 h r after the burst which had brightened by 
nearly a factor two 4.7hr later (ILevan et al.l 120071 ). The NIR afterglow was consistent with 
the location of the optical SDSS detection. Optical follow up de tected an object cons istent 
in location and brightness with the object seen in the SDSS data (IMalesani et al.ll2007l ). The 
NIR follow up was continued the following night when we triggered Target-of-Opportunity 
(ToO) observations with VLT/ISAAC under program 078.D-0416(E) (PI: Vreeswijk). A 
FORS 300V optical spectrum was obtained (program ID 078.D-0416(B)) at the position of 
the NIR afterglow and optical SDSS detection in order to secure a redshift measurement. A 
second ISAAC epoch was obtained on the following night to further characterize the NIR 
afterglow. Unfortunately, continued follow up was not possible as all available time for this 
run had been exhausted. Finally, late time NIR observations to secure estimates of the host 
were obtained with the WHT/LIRIS in the end of April and a K band point in January 
2008 as part of a GRB host program at VLT/ISAAC (program ID 177.A-0591(Q), PI: J. 
Hjorth). 

Using our best optical R band image and 30 objects from the USNO B1.0 catalog we 
determine the location of the host galaxy to be (a, 5) = (09 /t 52 m 23.30 s ±0'/20, +10°28 , 55"5± 
0'/12). Images of the GRB afterglow and host in R and K band are shown in Figs. [2] & [3j 
The afterglow in the K band has the same location within the errors, although a small offset 
of not more than a pixel (0'.'15) can be seen in Fig. [3l The host galaxy appears point-like in 
most of the filters, with the exception of the K band, where the deeper image reveals some 
faint extended emission. 



2.2.1. Imaging 



The optical data were processed in the standard way using IRAF@ and packages therein. 
The near-infrared data were processed using a combination of observatory pipelines (Eclipse, 
THELI) and IDL. Absolute photometric calibration was performed on t he best available 
imag ing data (R on May 7, 1 on Mar 12 and NIR on Mar 9) using SExtractor (IBertin fc Arnouts 
19961 ) and a fixed aperture of 3'.'0 radius. The zero point was established from standard star 
observations for the NIR (ISAAC) data, while for the FORS and NOT (optical) data we 
used six field stars from the SP SS data to photom etrically calibrate our R and I bands 
using the color transformations of iJester et al.l (120051 ). The remaining observations were cal- 



2 Image Reduction and Analysis Facility (IRAF), a software system distributed by the National Optical 
Astronomy Observatories (NO AO). 
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ibrated by offsetting the photometry relative to these reference images. The results are given 
in Table [U, and the reference images are denoted by daggers. The f oregro und reddening is 



estimated to E(B — V) = 0.0275 using the maps of ISchlegel et al.l (119981 ). but the column 
containing afterglow (Vega) magnitudes given in Table [T] have not been corrected for the 
Galactic extinction. 

Optimal sensitivity to variability is achieved by degrading each frame to the worst PSF 
quality in a series of images (for a given filter). The difference images are then inspected 
to identify residual emi ssion at the afterglow position. We performed image subtraction 
using the ISIS package ( Alard &: Lupton 19981 ) on all relevant imaging data. The reference 



frame for each filter (for registration and flux variability) was generally chosen as the latest 
epoch, but for J and H band these frames had a very low signal-to-noise and produced very 
poor image subtracted results. For these two series an earlier epoch had to be selected as a 
reference, although this was not ideal due to a possible afterglow contribution. The actual 
reference frames are indicated in Table [TJ Prior to running ISIS, the images were registered 
and transformed using the IRAF immatch package in order to ensure a common coordinate 
system and smooth processing by ISIS. We also copied an isolated star (a pseudo-PSF) of 
similar flux as that of the afterglow to an empty region of sky near the location of the GRB 
so that it would remain in the field after the image subtraction (for comparison purposes). 
In those frames for which we detected a residual (2a level) at the location of the GRB we 
employed fixed aperture photometry (using SExtractor) and an aperture radius of 1'.'5 to 
measure the relative flux between the residual and the reference star. In the opposite case, 
in which no residual (at the 2a level) could be detected, we estimated the detection limit 
by scaling our chosen PSF star by a decreasing amount and placing it at the location of the 
GRB afterglow/host. Image subtraction was then performed as described earlier and the 
magnitude difference was computed relative to the reference star (pseudo-PSF). This process 
was repeated until it could no longer be detected, determining the afterglow detection limit 
(or upper limit). 

The multi-color light curve (J, H, K bands) is shown in Fig. [1] and is characterized by 
an initial rise to peak during the first day in the K band and possibly also in the H band. 
The subsequent decline can be approximated by a power law, F oc t~ a . During the second 
and third days we fit the p ower law decline in t he K band with a = 2.1 ±0.4 typical for 



post-break lightcurves (e.g. lAndersen et al.ll2000l . their Fig. 4). No slope could be measured 



for the J and H bands. In the H band the afterglow is clearly detected in the first two 
epochs, but the poorer image quality and lower signal-to-noise of the last epoch prevented 
us from setting meaningful limits on the earlier epochs and we were therefore forced to use 
the third epoch as a reference image. In addition, the images originate from two different 
instruments, which clearly does not help. Consequently, we are unable to detect an afterglow 
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component at this epoch and may possibly overestimate the flux of the host (if there is a 
detectable contribution from the afterglow in this third epoch). We had similar problems 
with the R, I and J series, although the afterglow is not detected in these bands. In the J 
band there may be a marginal signal in the first epoch (at t + 0.19 days), but this is not 
significant enough to be measured. 

The K band light curve displays a prominent brightening duri ng the first day and 



although this has been observed in a f ew cases, this is a rare event (IPedersen et al.l 11998 



Stanek et al.ll2007l ; iGuidorzi et al.ll2007l ). There have been a number of explanations for such 



Tarn et al. 


2005; 


Lazzati et al. 


2002b) 



2.2.2. Spectroscopy 



On 2007 March 8 (starting at 02:38 UTC) using the VLT UT1/FORS2 instrument with 
the 300V grism we obtained three spectra with a total exposure time of 5400 s. The spectrum 
displayed a featureless continuum from 3700 to 9500 A, apart from a single emission line at 
~9306 A. This was attributed to [On] A3 72 7 at a redshift of 1.5 on account that no other 
emission or absorption lines were detected (IJaunsen et al.l 120071 ). The line significance is 
relatively high, but it is blended with the telluric skylines at 9306 A and onwards and with 
a resolution of about 12 A this prevented the line from being unambiguously identified. We 
have considered other potential line identifications such as H/5, [Om] A5007 or Ha, all of 
which are rejected on account that we expect to see additional and roughly similar bright 
lines in the wavelength range 6900-7140 A. 

In an attempt to resolve the [O n] doublet (expected line separation is 7 A) , we observed 
the host galaxy with the 1028z grism as part of our large program 177.A-0491(J) (PI: Hjorth). 
The observations were carried out in service mode on 2007 May 7 using VLT/FORS2 with 
a total exposure time of 2600 s. The seeing was Cf!7 and sky transparency was clear. The 
spectra were processed with special emphasis on obtaining a good sky subtraction. The issue 
is the fact that the skylines in these spectra are somewhat tilted, i.e. not perfectly aligned 
with the y-axis. In a traditional reduction the frames are first wavelength calibrated, which 
involves an interpolation of the original data values, and then sky subtracted. Since the steep 
flanks of the skylines are not sufficiently well sampled, the detailed shapes of the skylines 
in the interpolated frames vary from row to row, leading to systematic r esiduals i n the sky 
subtracted fram es. We dealt with thi s prob lem by using the method of iKelsonl (120031 ). as 
implemented in iMilvang- Jensen et al.l (120081 ). Since we are only interested in the limited 
region around the emission line, we manually inspected and cleaned cosmic rays only around 
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the line itself. A cosmic ray on the peak of the emission line of the first spectrum was cleaned 
by interpolating over the neighboring pixels. A new wavelength solution was obtained by 
identifying the numerous sky lines t hroughout the sp ectrum. The identification of sky lines 



was done using the UVES sky atlas (iHanuschik! 120031 ) . Finally, the new wavelength solution 



was applied to the stacked sky subtracted spectrum. 

The line spread function (LSF) is estimated from the sky lines, giving FWHM ~ 2.5 A, 
while the detected emission covers approximately 20 A from 9295 to 9315 A. If the emission 
line is due to the [O n] doublet then a peak to peak separation of nearly 7 A is expected 
at this redshift. We fit a Gaussian function to the blue component of the doublet with a 
fixed width corresponding to the LSF. The red component, which is severely affected by the 
skyline residuals, is assumed to be 30% brighter than the blue component and separated by 
7 A. Fig. H] shows the spectrum of the doublet with the profiles representing the unperturbed 
emission lines. The red component of the [On] doublet (vacuum rest-frame wavelengths 
A3727.09, 3729.88 A) is observed at wavelength 9302.6 ±0.2A, which corresponds to a redshift 
of z = 1.49594 ±0.00006. 

The flux calibrated 300V spectrum was normalized to the R band photometry. The 
integrated flux of the emission line was measured to be (1.3±0.2) x 10 -16 erg s _1 cm -2 . As this 
emission is significantly affected by the sky line residuals, we use our simple model to estimate 
a correction factor of 1.3 (by taking the ratio of the model and the spectrum). This is used 
to estimate the true flux of the [On] doublet, which becomes /[oii],obs = (1-7 ± 0.6) x 10~ 16 
erg s _1 cm -2 . Note that due to the uncertainty in the flux correction we have added a 
systematic error. T he star formation rate (SFR) can be computed using the relation to [On] 



luminosity given by iKennicuttl (119981 . Eq. 3), 

SFR [0 n] = (1.4 ±0.4) x 10- 41 L [O ii] M^r" 1 = 34+ 2 * M^r" 1 
where L[oii] = 47r<i 2 /[oii],obs and d\ is the luminosity distance. 



3. The host galaxy 

With an almost featureless optical spectrum, the host galaxy at first seems difficult 
to model. In the optical we have no indication of an afterglow contribution (variability) 
and the magnitudes are therefore taken as measurements of the host. The late time NIR 
photometry is also interpreted as a measure of the NIR host properties. This gives a total of 
10 independent magnitudes; SDSS /ugriz, FORS/i?, NOT// and LIRIS/ JHK. The adopted 
host magnitudes are listed in Table [2j 

The optical colors are consistent with the featureless spectrum, in that they seem to 



- 9- 



represent a nearly flat continuum. This trend does not continue at the very red end of 
the spectrum and in the NIR bands. The SDSS z band deviates from the other bluer 
bands, although with a large error, and the emission in the NIR peaks in the H band. We 
investigate the SED properties by exploring the possible significance of typical emission lines 
such as H a, H f3, [O n] and [O in] . The emission-line strengths depend primarily on star burst 
activity (star formation rate), the age of the stella r population, extinctio n and metallicity. 
The host SED is modeled using the ne wfl Hvper-Z faolzonella et aJ boOoh (vll) along with 
four spectral synth esis models based on I Coleman et al.l (119801 hereafter CWW) and two 
empirical templates iKinney et al.l (119961 ). The Hyper-Z code convolves the templates with 
the filter transmission curves to obtain colors which are compared to those observed by 
means of a minimum chi-squared. This provides some constraint on the host extinction, 
qualitative information on the SED and an independent constraint on the redshift (which 
for this reason is left as a free para meter). The CWW templates lack detailed emission 
features, so the IKinney et al.l (119961 ) star burst templates (SB1, SB2) are crucial in order 
to check thei r relev ance in explaining the observed SED. The difference between the two 
Kinney et al.l (119961 ) templates is primarily the amount of internal extinction. For the SB1 
template E(B - V) < 0.1 and for SB2 0.11 < E(B - V) < 0.21. A sati sfactory fit was 
achieved for all extinction curves we explored, but the MW (jSeatonl 1 19791 . hereafter S79) 
curve gave the best result with the IKinney et al.l (119961 ) SB2 template and a relatively well 



constrained redshift z = 1.55_ ; 26 and a \ /dof = 2.3/9. The chi-squared values were all 
admittedly low, probably suggestive of the errors being color large (possibly due to the SDSS 
errors) . 

From the spectroscopic results we have a firm constraint on the emission line and this 
can be implemented in the SED fitting as a pseudo narrow-band measurement at the wave- 
length of the detected emission line. The NB 930 7 filter has an effective wavelength of 9307.5 
A and a 19 A bandpass. The flux of the NB 930 7 filter was determined from the flux calibrated 
spectrum, giving NB 930 7 = 21.1 ± 0.5 mag (AB). Running Hyper-Z on the 11 magnitudes 
resulted in a best fit redshift of z = 1.50 ± 0.01 and x 2 /dof = 2.5/10 and the same combi- 
nation of the MW (S79) extinction curve and SB2 template. The resulting SED is shown 
in Fig. We thus conclude that the redshift of the host of GRB 070306 is z = 1.4959, as 
deduced from spectroscopy. Results from the host SED fits are summarized in Table [3J 

The broad band photometry fits to the SED templates do not indicate significant red- 
dening. For the various extinction laws we find the upper limit of Ay in the range 0.1 0.5. 

For the best fitting (SB2) template we find an upper limit of Ay < 0.45. 



http : //www. ast . obs-mip . f r/users/roser/hyperz/ 
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At the host redshift the [On] emission occurs in the range A 9304-9310 A, near the 
edge of the 1028z and 300V grisms. The emission line of [Nem] is expected at 9660 A at 
this redshift, but unfortunately falls outside the range of both grisms. The best option to 
independently verify the host redshift is therefore to obtain low resolution spectroscopy in 
J and H bands to search for the [O in] , H /3 and H a emission lines. 

Our best fittin g SEP template is a sta r burst galaxy template. This is typical among 
GRB host galaxies (IChristensen et al.l 120041 ) . The UV continuum flux in a star burst galaxy 
is dominated by the young (bursting) star population, and hence an appropriate measure of 
the star formation rate. We compute the SFR from the rest frame UV flux at 2800 A. The 
UV flux density is measured at the observed 2800(1 + z) A wavelength in the flux calibrated 
spectrum and gives /2soo,obs = 3.0±0.2 /xJv , which correspon ds to L v = 47r<i 2 /'2800,obs/(l+<2 ; ) = 
(17.1 ± 1.2) x 10 28 erg s" 1 Hz" 1 . Using iKennicuttl jl99fiL Eq. 1), which relates the SFR 
to the luminosity at A2800, we obtain SFR 2 8oo = 24 ± 2 M yr _1 . The error on the 
SFR does not include the roughly 30% error of the Kennicutt relation. The normalized 
luminosity L/L* is estimated from the rest frame absolute B magnitude, Mb = —22.3, 
which is computed by Hyper-Z. Adopting M* = —21.0 as the characteristic luminosity of 
field galaxies we find a normalized luminosity for the host of L/L* = 3.3. This is bright 
in co mparison to other G RB hosts, which primarily are sub-luminous (L < L*) galaxies 



e.g. 



Fruchter et al.l 120061 ) . The luminosity weighted star formation rate sSFR = 7.3 M 



yr 1 { L/L*) 1 . A similar U V estimate at a rest frame wavelength of 1500 A can be computed 
using iMadau et al.l (119981 ) and the observed flux density /i5oo,obs = 1-8 ± 0.6 /iJy. This gives 
L v = (10.3 ± 3.5) x 10 28 erg s" 1 Hz" 1 and an SFR 1500 = 14.4 ± 4.9 M yr" 1 . These SFR 
estimates are consistent with each other (bearing in mind the 30% uncertainty of the SFR 
relation) and with the SFR([Oirj) computed in Sec. 12.2.21 However, the SFR va lues are 
higher than most of the GRB host galaxies studied in e.g. IChristensen et al.l (120041 ) . 



4. Afterglow and extinction 

Variability in GRB 070306 was established early on thanks to the WHT/LIRIS data, 



(Rol et al. 


2007a: 


Levan et al. 


2007) 



ing two nights confirmed the K band variability, and revealed variability in the H band. The 
complete JHK light curve is shown in Fig. [U The variability of the afterglow component 
is most dramatic in the K band and less pronounced in the H band. In the J band only 
marginal variability, if any at all, is visible during the first night. An afterglow was thus 
detected in the H and K bands during the first three nights. 
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As described in Section 12.2.11 we have used image subtraction to measure the afterglow 
flux, and if undetected set upper limits to its contribution. To establish the SED of the 
afterglow, we computed the semi- simultaneous flux of the afterglow based on the observations 
carried out on the second night between 2007 Mar 08.054 UT and 08.070 UT. The ISAAC 
data in this period are superior in quality and therefore best suited for a contemporaneous 
SED analysis in addition to the FORS R band observation at Mar 08.104 to characterize 
the SED. We used the AB-offsets calculated by Hyper-Z for the R,I,J,H and K filters and 
found the values 0.23, 0.45, 0.94, 1.41 and 1.87, respectively. 

With an afterglow detection in only two bands and two additional upper limits there 
is a clear degeneracy in estimating the extinction and spectral slope. We need to introduce 
constraints on the spectral index and one option is to assume that the optical spectral index 
is identical to the X-ray, i.e. (3o = fix- We fit a number of extinction curves (see Tabled]) to 
the broad band optical/NIR data and compare the fits in terms of chi-squared (see Table H]). 
The fit is made by assuming an intrinsic (unobscured) spectral slope and fitting the amount 
of reddening (Ay) required such that the corrected points match the assumed spectral slope. 
The resulting Ay values are very similar and mostly consistent with each other within the 
errors. It is difficult to discriminate between the various extinction curves and the MW and 



LMC curves give the same results, but the SMC extinction curve of iPeil ( 119921 ) (Ry = 2.93) 
does provide a slightly better fit than the others. By requiring that the NIR slope matches 
the hardest X-ray slope (fio = 1-1), we obtain an Ay = 4.9 ± 0.6. However, the extinction 
corrected H and K points do not match the extrapolation of the X-ray power-law and a 
single power-law model is therefore excluded by this fit. If the errors are increased by a 
factor of three then a marginally consistent result can be ach i eved. The only extinction 



curve which does allow for the PL model to succeed is the IPeil (119921 ) SMC type curve. In 
the context of the synchrotron model, there is an additional possibility for the optical/X-ray 
spectrum, namely a broken power law with the cooling break between the two bands. This 
requires that fio = fix — 0.5 and gives a spectral slope of fio = 0.7 ± 0.1 using the fix value 
from Sec. 12.11 The best fitting extinction curve is obtained using Ay = 5.4 ±0.6. The errors 
represent the uncertainty related to the observational errors, which are dominated by the 
uncertainty of the H — K color. The SED and extinction curve fits are shown in Fig. [6J 

The constraints on the total extinction, A\, by exploring the possible X-ray spectral 
slopes and cooling break locations under the assumption that one of the two power-law 
models is correct. In the case of a single power-law we require that the extinction is not 
less than the shallowest X-ray slope supported by the data, ie.fi > 1.1 . Similarly, for the 
broken power-law model we use the maximum possible cooling break frequency to determine 
the minimum total extinction. The maximum total extinction comes from the limit of the 
minimum cooling break frequency, which is set by the frequency of the bluest band with an 
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optical/NIR afterglow detection (eg. the H band). The results are given in Tableland are 
visually indicated by the shaded regions in Fig. [6j 



5. Summary and discussion 



The observational data of GRB 070306 presented here suggests that this is a dark burst 
with a ve ry red NIR aftergl ow and that the host (z = 1.4959) is brighter than typical host 
galaxies (IBerger et al.l 120071 ). GRB 070306 possesses a higher than normal star formation 
activity in comparison to other GRB hosts. The afterglow SED analysis shows that signif- 
icant reddening is required to explain the observed optical/NIR afterglow assuming typical 
power-laws between the optical/NIR and X-rays. From the combined constraints of the X- 
ray and optical/NIR data and the assumption of a single or broken power-law, we can set 
limits on the total extinction without requiring a specific extinction law. We find that a 
single power-law is not supported by the data f or most of the extinction laws explored, while 
a broken power-law and Ay = 5.4 ±0.6 (SMC, |Peilll992l ) best describes the observed colors. 



The Galactic gas-to-dust ratio (IPredehl & Schmittl Il995l ) corresponds to N H = 1.7 x 



10 21 cm -2 /v4y (calibrated using X-ray absorption). Using the same ratio, the observed X-ray 
Nh would correspond to Ay = 15, which is much larger than observed, indicating a lower 
dust-to-gas ratio. Similar results have been found in the past, attribu ting the apparent deficit 



to dust destruction , dust depletion, or an unusual ISM composition (IGalama fc Wijersll2001 
Watson et allbood : ICampana et~al1l2007h . 



In the case of GRB 070306, the combination of what seems to be a highly obscured 
(and reddened) afterglow and a host with little or no indication of reddening illustrates the 
general concept that GRB progenitors reside in particularly dense regions within their host 
galaxies. If, by chance, the GRB line-of-sight traverses large columns of dust and gas while 
the host itself appears unobscured, this suggests the dust is unevenly distributed and likely 
concentrated in smaller regions or clouds within the host galaxy. Such patchy patterns are 
in effect also seen in our own Galaxy, alth ough it does not necessarily imply that dust and 
gas trace each other. iPadoan et al.l (120061 ) find evidence supporting the latter and suggests 
spatial fluctuations of the dust-to-gas ratio, even on very small scales, exist in the MW. If 
the dust obscuration occurs in smaller regions with supposedly varying optical depths, this 
may also have implications on absorption line studies and the metallicities resulting from 
them. There are several element abundance stud ies in the literature indicating dust depletion 
(e.g. ISavaglio fc Falll I2004J ; IWatson et al.l 120061 ) . but few (if any) cases providing a robust 



comparison of emission and absorption line induced metallicities. Such studies may therefore 
offer some hope in providing a better understanding of the isolated and global properties of 
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these galaxies. 



have been re 


ported 


(Levari et al. 


I2OO6: 



in the past (IMasetti et al. 
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Klose et a . 


2000h and 


more 


Rol et al. 


2007bl: 


Castro-Tirade- et al. 


2007 


r 


ranvir et al. 


2008) 



NIR observations are crucial in detecting the afterglow of dark bursts, and GRB 070306 

illustrates the importance of following NIR afterglows with 8 10m class telescopes. With 

the increasing number of apparent dust obscured (very red) afterglows, all of which are clas- 
sified as dark, we may be nearing an understanding of the fairly large fraction of dark GRBs. 
The nature of gamma-ray bursts and their association with regions of intense star-formation, 
and accordingly an abundance of gas and dust, provides a natural dust interpretation for 
dark bursts within the context of the synchrotron model. Although there are alternative 
explanations, the dust-related interpretation seems intuitively inline with the observed red 
color, the range of redshifts and host propertie s. An alternative i s that the X-ray and op- 
tical afterglows have separate physical origins (jPerley et al.l 120081 ) . As the number of well 
observed dark bursts increase we are likely to discover many more of these systems, and 
hopefully settle the ambiguities of interpretation. 
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Table 1. Observing log and results from photometry 



Date 


Time a 


Instrument 


Band 


Seeing 


Tint 


Magnitude 13 


Error 




Error 


UT 


days 








S 


Vega 


mag 




uJy 


Mar 06.8889 


0.191 


WHT / LIRIS 


J 


1"1 


1200 


> 22.4 




< 1.7 




Mar 08.0542 


1.357 


VLT / ISAAC 


J 


0"9 


1040 


> 22.0 




< 2.5 




Mar 09.2342 


2.537 


VLT / ISAAC 


J 


077 


2560 


(20.90)t 


0.13 


6.82 


0.75 


Apr 30.8871 


55.190 


WHT / LIRIS 


J 


0"8 


1800 











Mar 06.8670 


0.170 


WHT / LIRIS 


H 


17 1 


1200 


19.78 


0.12 


12.31 


1.36 


Mar 08.0836 


1.386 


VLT / ISAAC 


H 


078 


1152 


20.52 


0.21 


6.23 


1.20 


Mar 09.2074 


2.510 


VLT / ISAAC 


H 


0"6 


1344 


(20.00)+ 


0.20 


10.05 


1.48 


Apr 30.9242 


55.227 


WHT / LIRIS 


H 


l'.'O 


1800 










Mar 06.8416 


0.144 


WHT / LIRIS 


K 


l'.'O 


1500 


18.28 


0.15 


31.92 


4.41 


Mar 07.0361 


0.339 


WHT / LIRIS 


K 


079 


1200 


17.61 


0.14 


59.16 


7.62 


Mar 08.0701 


1.373 


VLT / ISAAC 


K 


0"7 


1152 


18.53 


0.15 


25.35 


3.50 


Mar 08.2153 


1.512 


VLT / ISAAC 


K 


0'/6 


1152 


18.50 


0.15 


26.06 


3.60 


Mar 09.1879 


2.490 


VLT / ISAAC 


K 


076 


1152 


19.83 


0.19 


7.66 


1.34 


Apr 30.9699 


55.272 


WHT / LIRIS 


K 


079 


1800 










Jan 30.2244 


329.526 


VLT / ISAAC 


K 


075 


1920 


(19.60)t 


0.19 






Mar 06.9523 


0.255 


NOT / StanCam 


I 


l'.'O 


1200 










Mar 12.0615 


5.364 


NOT / ALFOSC 


I 


076 


3600 


(22.15) 1 " 


0.19 


3.57 


0.62 


Mar 08.1042 


1.407 


VLT / FORS2 


R 


079 


120 


> 23.5 




> 1.3 




Apr 20.0507 


44.364 


VLT / FORS2 


R 


079 


180 


> 23.5 




< 1.3 




May 07.0040 


61.306 


VLT / FORS2 


R 


077 


60 


(22.78)+ 


0.23 


2.53 


0.54 



a Time since burst trigger t = 2007 Mar 06.6976 UT 

b The photometry was obtained by measuring the flux in circular apertures of 3'/0 radius (using SExtractor) on the 
best data available (around Mar 09.2, t=1.38 days and best in terms of sub-arcsecond seeing, photometric conditions and 
telescope aperture). All other magnitudes and upper limits are measured from image subtracted frames and calibrated 
using the reference frames, which are identified by the daggers (see Sec. 12.2.11 for details). 

°The flux values have been co mputed using the AB offsets (see Section H]) and correcting for Galactic extinction 
according to iSchlegel et al.l (|1998l ) 
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Table 2. Host magnitudes 





u 


g 


r 


i 


z 


R 


I 


NB9307 


J 


H 


K 


Mag(AB) 


23.17 


23.06 


22.96 


22.84 


21.90 


22.98 


22.62 


21.1 


21.84 


21.41 


21.47 


Error 


0.46 


0.17 


0.25 


0.40 


0.54 


0.23 


0.19 


0.5 


0.12 


0.16 


0.20 



Note. — The SDSS ugriz colors are provided by lCool et al.1 (|2007l ). See Section[3]for a description 
of the pseudo narrow-band NB9307. 



Table 3. Results from the GRB 070306 host galaxy SED fit of ten broad-band filters and 



one pseudo narrow-band filter using Hyper-Z. T 
to those implemente d in Hyper-Z, in addition to 



(SMC); lAUenl (ll976T) f or the Milky Way (MW), ISeatonl (1 19 791 ) fit bv iFitzpatrickl (119861) fo r 



le e xtinct ion curves used here correspond 
Pea dl992h fo r the S mall Magell anic Cloud 



the MW. iFitzpatrick (119861) for the Large M agellanic Clou d (LM C), IPrevot et al.l (119841 ) 
and iBouchet et al.l (119851 ) for the SMC and ICalzetti et al.l (120001 ) for starburst galaxies 

(SB). 



Ext. law 


SED 


z 




A v 


max(Av) a 


^low 




MW (A76) 


SB2 


1.49 


0.256 


0.17 


0.43 


1.49 


1.51 


MW (S79, F86) 


SB2 


1.50 


0.249 


0.17 


0.45 


1.49 


1.51 


LMC (F86) 


SB2 


1.49 


0.268 


0.17 


0.45 


1.48 


1.51 


SMC (P84, B85) 


SB1 


1.49 


0.316 


0.03 


0.23 


1.48 


1.51 


SB (COO) 


SB2 


1.49 


0.269 


0.23 


0.47 


1.48 


1.51 


SMC (P92) 


SB1 


1.49 


0.325 


0.00 


0.08 


1.48 


1.51 



^The (la) upper limit of Ay. 
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Table 4. Rest-frame Ay and x 2 (dof = 1) results from fitting the afterglow SED for 
various extinction curves. While the MW and LMC curve s provided similar results, the 
overall best fitting curve was the SMC curve of iPeil (119921 ). The da ta are not consisten t 
with th e (single) PL model, e xcept in the case of the SMC curve of iPrevot et al.l (11984 ): 



Bouchet et al.l (119851 ). where a marginal consistency could be achieved. 



Power-law Broken power-law 



Ext. law 



Ay X 2 M 



X 



MW (A76) 


4.81 


0.17 


5.31 


0.27 


MW (S79,F86) 


4.81 


0.17 


5.31 


0.27 


LMC (F86) 


4.81 


0.17 


5.31 


0.27 


SMC (P84, B85) 


6.00 


0.65 


6.71 


1.13 


SB (COO) 


5.15 


0.48 


5.75 


0.81 


SMC (P92) 


4.93 


0.17 


5.45 


0.26 



Table 5. Total extinction constraints. 



Power-law Broken power-law 



Filter mm(A\) max(A\) min(A\) max(A\) 



R 6.0 ... 3.8 

J 6.0 ... 4.4 

H 5.4 6.9 2.8 6.9 

K 4.2 5.7 1.5 5.6 
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1000.0 



'(f) 
c 

CD 

X 
3 



100.0 
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Fig. 1. — The multi-color lightcurve (J, H, K) of GRB 070306 and estimated late-time decay 
index, a — 2.1 ± 0.4. The vertical dashed line marks the simultaneity time, t + 1.38 days. 
Upper limits are denoted by an arrow. 
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OA+host host OA R 

Fig. 2. — The R band image at Mar 08.10 (left), combined Apr 20 and May 07 (middle) 
and the difference image (right) obtained using ISIS. No residual emission is detected and 
the estimated upper limit of the afterglow is R > 23.5 (estimated as described in Sec. I2.2.1T) . 
North is up, and East is left. 




OA+host host OA K 

Fig. 3. — Image subtraction (ISIS) is used to detect variability and in this montage we show 
the K band image of Mar 08.21 with contribution from the afterglow and host (left). In 
the middle frame the late-time (host) image from 2008 Jan 30.22 is shown, clearly revealing 
its extended nature. To the right the difference image reveals the significant contribution of 
the afterglow of GRB 070306 in the K band. The images, which have been registered using 
field stars, show that the afterglow is centered on the brightest region in the host galaxy 
(although a small offset towards the north-west may be present on the order of a pixel, 0'.'15). 
North is up, and East is left. 
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Fig. 4.— The [On] doublet A3726.1, 3728.8 A at z = 1.4959 of GRB 070306 observed on 
2007 May 7 using the VLT/FORS2 1028z-grism. The red component is affected by skyline 
residuals in the 9306 - 9324A wavelength range and this is evident when comparing the 
spectrum to the prominent features in error spectrum. A Gaussian function with the same 
width as the LSF was fit to the blue component and the red component was assumed to be 
30% brighter than the red and the peaks separated by 7 A. The result is shown by the solid 
lines and the combined profile by the dotted line. 



-27- 



100.0 F 



0.0 



.0 



0. 








R I NE 
u g r i z 



'9307 



Js 



Ks 



5.0x10 3 1.0x10 4 1.5x10 4 2.0x1 4 2.5x10 4 



wavelength (A) 



Fig. 5. — Best fitting SED (% 2 = 2.5/10) using Hyper-Z and a iKinney et al.l (119961) type 
SB2 template at redshift z = 1.50 with a MW-type extinction curve ( Seatonll979 ) and 
Ay = 0.17. The fit is based on 10 optical and NIR broad-band filters and an additional 
pseudo narrow-band filter NB 9 3 07 . 
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Fig. 6. — Various extinction curves are fit to the optical and NIR afterglow measurements 
at to + 1-38 days by assuming synchrotron emission in the form of a single power-law (PL) 
or broken power-law (BPL) spectral slope. A separate power-law fit to the XRT data gave 
a spectral slope of (3x = 1.2 ± 0.1, which is used to constr ain t he op tical/NIR slope, flo- 
The best fit is obtained using the SMC extinction-curve of IPeil ( 19921 ) (represented by the 
restframe visual extinction, Ay). The observed optical/NIR measurements are shown as 
large filled circles, whereas the large open circles represent the same data corrected for the 
best fit extinction. The small filled dots with error bars represent the XRT data corrected 
for a fixed Galactic absorption, while the small open dots represent the same data corrected 
for a variable equivalent hydrogen column density. The shaded areas represent the limits on 
the total extinction without assuming an extinctio n curve (s ee text). The single power-law 
model (left) is not c onsistent with the d a ta for this |Pe i (1992) SMC curve, but is marginally 
consistent using the Prevot et al. (jl984 ); Bouchet et al. ( 1985 ) SMC curve. 



